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1.0  Introduction 

The  structure  of  turbulence  within  and  near  inversion  layers  as  well  as  the  motion 
within  the  interior  of  an  outer  stratified  layer  are  of  interest  in  understanding  the 
communication  between  stratified  and  turbulent  layers.  Such  an  understanding  is 
imperative  in  predicting  the  penetration  of  a  turbulent  layer  into  a  contiguous  stratified 
region  (entrainment  phenomena)  and  the  dispersion  of  species  in  the  turbulent  layer.  A 
stunning  example  in  this  context  is  the  fluid  motion  in  the  lower  atmosphere  generated  by 
convective  turbulence.  Here  the  stratification  consists  of  three  layers;  the  bottom  layer  near 
the  surface  is  well  mixed  by  convective  turbulence  and  the  outer  weakly  stratified  layer  is 
separated  from  the  bottom  layer  by  a  strong,  stable,  somewhat  thin  layer  called  the 
inversion  layer.  The  height  of  the  bottom  layer  is  referred  to  as  the  mixing  height.  Thus 
the  buoyancy  frequency  profile  in  the  lower  atmosphere  can  be  approximated  as 

N(z)  =  0  for  0  <z  <  hj 

=  Ni  for  hj  <  z  <  hj  +  h 
=  Nu  for  hi  h  <  z  <  h]  +  h  +  h2 

where  z  is  the  vertical  coordinate  measured  from  the  ground,  hj  the  mixing  height,  the 
buoyancy  frequency  within  the  inversion  layer  of  thickness  h,  and  the  buoyancy 

frequency  above  the  inversion.  It  has  been  demonstrated  that  mean  wind  effects  such  as 
shear  production  of  turbulence  can  be  neglected  except  very  near  the  ground  (Deardroff 
1980)  and  hence  the  forcing  of  the  inversion  can  be  treated  as  due  to  shear-free  turbulence. 
This  is  somewhat  different  from  the  oceanic  case,  where,  albeit  a  three  layer  stratification 
consisting  of  an  upper  mixed  layer,  a  thermocline  and  a  deep  weakly  stratified  layer  can  be 
identified,  the  turbulent  production  at  the  base  of  the  mixed  layer  is  significant. 

The  work  reported  here  deals  with  the  interaction  of  turbulence  and  stratification  in 
a  three  layer  system  similar  to  that  described  above  and  the  nature  of  shear-free  turbulence 
near  sharp  interfaces.  Summaries  of  the  two  projects  are  given  below. 
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2.0  Experiments  with  Three<Layer  Stratified  Fluids 

The  experiments  were  carried  in  a  traditional  mixing  box  similar  to  that  used  by 
Turner  (1968).  For  experimental  convenience,  the  upper  layer  was  maintained  turbulent 
and  the  layers  below  were  stably  stratified.  Shear  free  turbulence  induced  by  mechanical 
means  (an  oscillating  grid)  was  used  to  drive  the  upper  turbulent  layer.  A  schematic  of  the 
experimental  apparatus  is  shown  in  Figure  1(a)  and  the  stratification  used  is  shown  in 
Figure  1(b).  Based  on  linear  internal  wave  theory,  one  expects  that  waves  generated  at  the 
inversion  layer  will  have  frequencies  satisfying  O)  <  Ni  and  the  waves  radiated  into  the 
outer  layer  will  have  frequencies  co<  .  Since  Ni>N^,  one  further  expects  that  waves  in 
the  inversion  satisfying  N^<  (o<  will  be  trapped  inside  that  layer  and,  owing  to  the 
build  up  of  energy,  break  thus  dissipating  energy.  During  the  wave  breaking,  turbulence  is 
generated  which  leads  to  localized  mixing  within  the  inversion,  reduction  of  the  buoyancy 
gradient  within  the  interface  and  growth  of  the  mixed  layer. 

A  three  layer  stratification  was  set  up  in  the  experimental  apparatus  using 
conventional  techniques.  The  laser-induced  fluorescence  (LEF)  technique  and  conductivity 
measurements  were  used  to  examine  the  wave  activity  within  the  stratified  layers.  A  two- 
dimensional  image  processing  technique  was  developed  to  probe  intcrfacial  mixing  events. 
Although  a  laser-based  local  Richardson  number  probe  was  used  to  detect  waves  in  the 
stratified  layers,  measurement  of  wave  parameters  using  this  method  was  difficult  owing 
to  the  low  wave  activity  in  the  stratified  layers. 

The  LIF  flow  visualization  revealed  the  existence  of  different  mixing  mechanisms  in 

different  Richardson  number  (Ri)  ranges.  Here  the  Richardson  number  is  defined  as  Ri  = 
,  where  and  Uff  are  the  length  and  velocity  scales  of  the  (undistorted) 

turbulence,  respectively.  At  low  Ri,  say  Ri  <  25.  the  interfacial  mixing  takes  place  owing 
to  the  scouring  of  the  interface  by  the  energetic  turbulent  eddies  in  the  mixed  layer.  As  Ri 
increases,  wave  breaking  becomes  important  and,  at  Ri  >  40,  mixing  is  dominated  by  the 
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breaking  of  internal  waves.  At  very  high  Ri,  the  interface  was  found  to  be  rather  calm  and 
non  entraining;  here  the  interface  is  dominated  by  molecular  diffusive  processes. 

The  maxima  of  the  internal  wave  energy  spectra  lie  between  the  buoyancy 
frequencies  of  the  inversion  layer  Ni  and  the  weakly  stratified  layer  Nu\  resonant  trapped 

waves  within  the  inversion  could  be  identified.  Measurements  of  these  wave  frequencies 
showed  good  agreement  with  a  theory  due  to  Carruthers  and  Hunt  (1992).  The  interfacial 
layer  thickness  was  measured  using  two  techniques.  The  conductivity  measurements 
showed  that  the  normalized  interfacial  layer  thickness  is  a  slowly  varying  function  of  the 
Richardson  number,  i.e.,  dlL^j  -  Ri  The  concentration  measurements,  on  the  other 
hand,  showed  a  much  more  rapid  decrease  with  Ri;  i.e.,  dlL^  ~  Ri  The  r.m.s. 
amplitude  of  the  normalized  interfacial-wave  distortions  was  found  to  decrease  with  Ri  as 
\C  I  iL^j  ~  Ri  and  the  energy  spectrum  of  the  interfacial  distortions  showed  the 
existence  of  a  region  with  a  -2  slope,  in  accordance  with  the  predictions  of  Moffatt  (1984) 
and  Fernando  &  Hunt  (1992).  The  r.m.s.  fluctuation  of  normalized  interfacial  vertical 
velocity  was  found  to  be  a  decreasing  function  of  Ri  as  (w^f  -  Ri 

Complete  description  of  the  experiment  and  the  experimental  results  will  be  given  in 
Perera,  Fernando  &  Boyer  (1992). 

3.0  Modelling  of  Turbulence  Near  Sharp  Interfaces 

A  model  for  turbulent  mixing  across  sharp  density  interfaces  subjected  to 
homogeneous  turbulence  was  developed  in  this  study.  Th''  analysis  was  essentially  based 
on  the  rapid  distortion  theory  approach  of  Hunt  (1984),  which  reckon  on  the  empirical 
observation  that  the  dissipation  rate  of  turbulent  kinetic  energy  above  an  interface  is 
independent  on  the  height  The  homogeneous  turbulence  was  modelled  using  the  classical 
von  Kaiman  spectra  and  the  interfacial  motions  were  assumed  to  be  governed  by  the  linear 
internal-wave  equations.  It  was  further  assumed  that  the  first  mode  of  internal  waves 
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governs  the  interfacial  motions,  thus  simplifying  the  interfacial  boundary  conditions;  the 
latter  has  been  a  major  bottleneck  in  modelling  fluid  motions  near  sharp  density  interfaces. 

Two  configurations,  namely,  (i)  a  turbulent  fluid  layer  separated  from  a  non* 
turbulent  heavy  fluid  layer  by  a  density  interface  (single-sided  stirring)  and  (ii)  a  density 
interface  sandwiched  between  two  layers  of  equal  turbulent  intensities  and  length  scales 
(double-sided  stirring),  were  considered.  The  former  case  is  applicable  to  situations  of 
mixed  layer  deepening  in  stably  stratified  fluid  masses  due  to  convective  stirring  or  wave 
breaking  whereas  the  latter  is  important  in  analyzing  the  interfacial  migrations  in  density 
step  structures.  Detailed  measurements  on  the  interfacial  structure  for  the  case  (i)  have  been 
made  by  Hannoun  &  List  (1988),  when  the  interface  is  forced  by  oscillating- grid  induced 
shear-free  turbulence.  These  results  were  effectively  utilized  to  obtain  a  physical  insight  for 
the  problem  and  for  the  validation  of  the  model 

The  existence  of  resonant  modes  at  the  interface  was  identified  when  co  <  Ri/2;  the 
waves  at  or  near  resonant  frequencies  amplify  during  the  evolution  and  break  down  owing 
to  local  shear  instabilities  to  form  isolated,  intermittent,  patches  at  the  interface  (here  Ri  is 
based  on  the  buoyancy  jump  at  the  interface).  The  patches  act  as  energy  sinks,  as  a  result 
of  mixing  and  dissipation,  and  are  responsible  for  the  maintenance  of  internal-wave  energy 
at  a  finite  level.  Such  breaking  events  cannot  be  described  by  linear  theory;  modelling 
based  on  physical  and  mathematical  reasoning  was  introduced  to  close  the  problem.  Based 
on  the  experimental  observation  that  the  time-traces  of  the  interfacial  displacement  at  a  given 
location  has  discontinuities,  the  following  spectral  form  was  proposed  for  the  vertical 
velocity: 


V33  (0),o)  = 


^co)  formxuk; 

l{oic)  foroxco^ 


where  I(o))  is  the  spectral  form  predicted  by  the  linear  theory  and  is  the  frequency  below 
which  non-linear  breaking  of  waves  becomes  important.  The  cut  off  frequency  0)^  was 
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evaluated  using  the  physical  argument  that  the  ratio  of  the  energy  G  contained  in  the  vertical 
velocity  spectra  in  non-linear  and  linear  regimes  is  constant;  this  assumption  led  to  the 
condition  m^  =  Ri/2ca^  is  a  constant,  independent  of  the  Richardson  number. 

It  was  found  that  the  vertical  velocity  data  give  excellent  agreement  with  the 
predictions  for  case  (i),  if  m^  »  0.5  (corresponding  to  G  =  1.25)  is  chosen.  The 
calculations  for  various  quantities  were  carried  out  using  m^  =  0.5  for  both  cases  (i)  and 
(ii).  In  calculating  the  displacement  spectra  using  the  modelled  spectra,  the  low  frequency 
end  near  O)  =  0  behaves  as  co'^,  and  thus  the  existence  of  a  low  frequency  cut  off,  dictated 
by  the  lowest  frequency  of  the  motion  in  the  mixed  layer,  was  suggested.  The  proposed 
form  of  the  displacement  spectra  was 

V*(to.O)  = 

Based  on  previous  work,  the  non-dimensional  (Uo  was  selected  as  0.18.  Using  the 
above  spectral  forms  with  nic  =  0.5,  coq  -  O.IS  and  the  standard  von  Kaiman  form  of 
spectra  for  homogeneous  turbulence,  the  Richardson  number  Ri  dependence  of  non- 

_ _ .j,2 

dimensional  r.m.s.  vertical  velocity  (w^)  ,  displacement  1C  )  and  pressure-velocity 
correlation  Apw  at  the  interface  were  calculated  as  follows: 


^33 (m.oyto^  form >0)0 
I V33 (o)o,0yco5  for o)< 0)0 


Single-Sided  Stirring 


Two-Sided  Stirring 


=  1.8  Ri-»  ” 
=  3.2  Ri-“  ” 


Apw  =  5.0  Ri 


-0.49 


=  3.0  Ri-O  ” 

=  5.4  Ri-O  ” 
Apw  =  16.7  Ri-®-^’ 


However,  it  was  pointed  out  that  the  form  of  joint-wave  number  frequency  spectra 
Xij(k,co)  of  homogeneous  turbulence  used  in  the  present  work,  which  is  an  essential  input 


5 


to  the  modelling,  does  not  yield  an  Eulerian  frequency  spectra,  consistent  with  the 
experimental  observations,  upon  integradon  over  the  wave  number  space.  Calculations 
performed  with  appropriate  modifications  showed  that  the  predictions  based  on  the 
assumption  G  =  1.0  (i.e.,  1x1^  =  0.36)  also  give  a  good  agreement  with  the  experimental 
data.  The  resulting  predictions  are  as  follows: 


Single-Sided  Stirring  Two-Sided  Stirring 


=  1.8  Ri-O-33 

=  3.2  Ri-O  ” 

=  3.2  Ri-O-83 

(3^^  =  4.9  Ri-o  ” 

Apw  =  5.0  Ri-o-^^  Apw  =  1 1.8  Ri-®-*’ 


An  entrainment  model  was  developed  to  calculate  entrainment  rates  across  density 
interfaces.  The  model  predictions  were  found  to  be  in  good  agreement  with  the  entrainment 
data  obtained  in  previous  oscillating-grid  experiments  such  as  Turner  (1968)  and  E  and 
Hopfinger  (1986).  Turner  (1968)  conjectured  that  the  entrainment  rates  for  single  and 
double-sided  cases  ought  to  be  the  same.  One  of  the  important  findings  of  the  present 
analysis  is  that  these  two  cases  are  different  and  the  different  entrainment  rates  observed  by 
Turner  (1968)  can  be  real.  The  differences  arise  from  the  fact  that  the  profiles  of  buoyancy 
flux  and  the  interfacial  velocities  within  the  interface  are  dissimilar  for  the  two  cases; 
however,  the  power  law  dependencies  of  various  quantities  on  the  Richardson  number  Ri 
were  found  to  be  similar. 

The  present  analysis  is  more  complete  than  the  previous  theories  on  turbulent 
mixing  across  interfaces.  Hannoun  and  List  (1988)  used  the  analysis  of  Phillips  (1977), 
based  on  first  mode  of  internal  waves  at  sharp  interfaces,  to  interpret  their  results  but  the 
latter  does  not  take  into  account  the  effects  of  turbulence  on  the  interface  which  lead  to 
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resonant  breaking  of  waves.  Mory  (1990)  assumed  that  Kelvin-Helmholtz  instabilities 
resulting  from  the  sloshing  motion  of  the  eddies  on  the  interface  leads  to  mixing.  An 
approximate  analysis  indicates  that  such  a  mechanism  is  unlikely  to  play  an  important  role 
in  mixing,  except  at  small  Ri. 
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Figure  1.  (a)  A  schematic  of  the  experimental  apparatus  (b)  the  three  layer  stratification 

usal  for  the  experiment  The  positions  of  the  conductivity  probes  are  given  by 
Cl  C2,  C3  and  C4  and  the  traversing  conductivity  probe  is  represented  by 

Cs- 


